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Abstract 



•rH . Following Bogomolnyi's classical treatment of vortices, we develop a method 

>< 

for finding rigorous lower bounds to the Landau-Ginzburg free energy describ- 
ing unconventional superconductors in the absence of external magnetic fields. 
This allows a more precise description of the magnetic instabilities previously 
considered in these systems. In particular, we derive new sufficient conditions 
for the stability of both the homogeneous and inhomogeneous equilibrium 
states. 
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In recent years, a number of authors have considered the possibility of magnetic insta- 
bilities in unconventional superconductors. Within the framework of a phenomenological 
Landau-Ginzburg description of these materials |]J, a multicomponent order parameter can 
give rise to a number of couplings which make the system potentially unstable against spatial 
modulations. The stability of homogeneneous phases was first considered by Zhitomirskii 
||. Soon after that, Palumbo, Muzikar and Sauls fU described a magnetic instability in the 
absence of externally applied magnetic fields, associated with the creation of supercurrents 
due to spatial distorsions of the order parameter. They found the instability by studying 
the change in the free energy under small perturbations about the homogeneous phase. 

Using similar techniques, Zhitomirskii 0] later identified several possible inhomogeneous 
phases; he also derived some necessary conditions for the Landau-Ginzburg free energy 
functional to be bounded from below. This is a crucial question because the existence of an 
unstable mode, while automatically implying the instability of the homogeneous phase, does 
not guarantee the existence of an inhomogeneous phase. It could very well happen that the 
free energy is unbounded from below. 

On the other hand, it is important to realize that sufficient conditions for the instability 
of a given phase are not the same as sufficient conditions for its stability. The standard 
analysis of small perturbations can only provide information about the former. The latter 
can sometimes be derived from a study of the minima of the potential provided the gradient 
terms are positive definite. Both methods tend to give conservative bounds, which often 
leave a region of phase space intractable. Our aim here is to bridge the gap between the two 
and obtain information on the stability in these intermediate regions. 

In this paper we present a general method which can be used to obtain new sufficient 
stability conditions for Landau-Ginzburg free energies describing unconventional supercon- 
ductors. The basic idea, inspired by Bogomolnyi's classical treatment of vortices |J, is to 
rewrite the free energy in such a way that lower bounds to its value are apparent. This is 
achieved by using the boundary conditions to transform and recombine different terms into 
perfect squares, so that positivity becomes obvious. As we shall see, our results are comple- 
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mentary to previous studies based on the analysis of small perturbations to the homogeneous 
equilibrium states, which yield necessary stability conditions. In some cases we are able to 
establish both necessary and sufficient conditions, thus providing the exact location of the 
phase transition. 

Following || and || we consider an order parameter V = (Vi, V2) transforming according 
to the two-dimensional E lg representation of the point group D sh , but our method easily 
generalizes to other cases. We will consider z-independent configurations, with B = d x A y — 
d y A x and D = V — %A. The free energy density can be written as 

F = F B + F G + F P (1) 

where Fb = \B 2 is the magnetic energy contribution (in the absence of externally applied 
field), Fq depends on the gradients, 

F G = KiDj/pDaVp + K 2 D a V a D p Vp + KzD a VpDpV a (2) 

and Fp is the potential 

Fp = ^P(VV-l) 2 + ^VxV) 2 (3) 

with VV = V1V1 + V2V2, and i(Vx V) = i(ViY2 — V2V1). Note that we have already performed 
several rescalings in order to bring the free energy into this particularly simple form. We 
must take (3 > and a > —f3 for global stability (these are necessary conditions). 

Let us first consider the homogeneous phases. Setting DV = B = we have to distin- 
guish two cases, depending on the value of a: 

• For a > 0, Fp is the sum of two positive squares and takes its minimum value Fp = 
for V real (iV xV = 0) with VV — 1, i.e., 

/cos#\ 

Vo = , o > (4) 

V sin# / 

This is sometimes referred to as the time reversal symmetric phase in the literature. 

• For —j3 < a < 0, (§) can be rewritten as 
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Fp = \ {0 + a) (yy _ _1_) 2 _ I ff _ (i p x VT) + i^L (5) 

and we have Fp > \-§t^- Fp attains its minimum value for \VV\ = \iV x V\ = ~Jt^, i-e., 



P 
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also known as the time reversal symmetry breaking state. Eqns. (f|) and (§) with A = 
represent two types of homogeneous phases. In order to study their stability, we must take 
into account the behaviour of Fq. Consider first configurations with A = 0. The contribution 
to Fq from a plane wave with V 2 = and k = (k x , k y ) will be proportional to 

(Kmfc' + KiA^lKl 2 (7) 

and this will be positive for k 123 = k% + k 2 + % > and «i > 0. These necessary conditions 
for global stability of F have been considered in || and 0, and in what follows they will 
be assumed to hold. 

Consider now general configurations with A 7^ 0. Fq can be written 

F G = Kiaafl Wi| 2 + I W 2 | 2 ) + ^(L^ViA^ + D^A^a) (8a) 

+ k 1 (\D 1 V 2 \ 2 + |iW| 2 ) + KaCA^DaVi + D x V 2 D 2 Vx) (8b) 

and it is clear that fl3a]) and ( [SbD will be positive definite forms in {DiVi, D 2 V 2 } and 
{-D1V2, A^i} respectively if 

|«2| < ^123 and |« 3 | < «i (9) 

If (^|) is satisfied we have > 0. Since the minimum value of Fp is obtained for homo- 
geneous configurations (with A = 0), it is clear that (|9]) represents sufficient conditions for 
the stability of the homogeneous phases. These conditions were given without proof in 0. 

We will now show that the domain of existence of the homogeneous phase is in fact larger 
than conditions (§) suggest. Consider the case a > 0. There are two key observations which 



allow us to enlarge the region of stability given by (|9]). The first one is that, up to integration 
by parts (which is allowed, for instance, by periodic boundary conditions) 

(D 1 V 1 D 2 V 2 + D 1 V 1 D 2 V 2 ) - (D 1 V 2 D 2 V 1 + D X V 2 D 2 V X ) =i(VxV)-B 

and we may write 

F g (ki,k 2 ,k z ) = F g {ki,k 2 - 5k,k z + 5k) + i5K(V xV) ■ B . (10) 

The second observation is that the extra term in (|H]) can be combined in the following way 
with Fp and F B : 

F B + F P + iSK(V xV)-B = ^[B + iSK(V x V)f 

+~(3(VV - l) 2 + ~(<r - 5k 2 MV x V)f . (11) 

Since this is positive definite as long as \Sk\ < yfa, this freedom in the choice of 5k in 
([10|) can be used to transform @ into 

|«2| < «i23 + \fv and |/c 3 | < k x + yfa (12) 

Basically we have "pushed" the boundary of the stability region given by (|9|) along vectors 
±(y / o : , —yfa). Within the region given by (Jl2|), F > 0, and the minimum value F = is 
attained by the homogeneous state (^). Thus the homogeneous phase is stable if ( |12] ) is 
satisfied 

We have just shown the sufficient character of fll2|) . Equations (|12|) are also necessary 
conditions. This follows from the consideration of perturbations V — Vq + 5^ to the homo- 
geneous state (U). Using rotational symmetry to set 9 = 0, one finds that F is unstable 
against perturbations of the form: 

( ° \ 

5V oc ; \k 2 \ > k 123 + y/a (13a) 

\ i cos ky J 



SV oc ( ] ; |k 3 | > Ki + v 7 ^ (13b) 

z cos kx 



These perturbations were considered by Zhitomirskii H, who showed that (12) are necessary 
conditions for the stability of (^). The fact that these are identical with our new sufficient 
conditions means that eqs. ( |T2]) define exactly the domain of existence of the homogeneous 
phase 

Now let us consider what happens outside this domain. The point is that finding an 
instability of the homogeneous phase does not in itself prove the existence of an inhomoge- 
neous phase. To do this, one also has to show that the free energy remains bounded below, 
and this cannot be established by studying its changes under small perturbations. 

Notice, however, that our method provides such a proof: according to @ Fp is bounded 
from below as long as a > —(3, therefore we may choose \8k\ < y/]3 + a, and conclude that 
inside the region defined by 

\k 2 \ < Km + \/f3 + a, \k 3 \ < Ki + \J (3 + a, (14) 

the total free energy density is bounded from below. This is obviously true irrespective of 
the sign of a, and in the case where a > it proves the existence of a stable inhomogeneous 
phase without resorting to numerical techniques. 

Note that eqs. ( |13aj ) and ( |13b| ) can be described as "longitudinal" and "transverse" 
modulations respectively and that for K3 < —Ki both perturbations are simultaneously 
unstable. In that region we may thus expect the occurrence of 2-D modulated solutions. 
This seems to be confirmed by the numerical simulations in ||. 

To summarize, we have suggested a new approach to the study of lower bounds on the free 
energy of condensed matter systems inspired on Bogomolnyi's ideas. Our basic observation 
is that integration by parts generally results in a mixing of the gradient, potential and 
magnetic contributions to the free energy which can be exploited to obtain sufficient (rather 
than necessary) conditions for stability. This gives complementary results to the standard 
analysis of small perturbations and we have shown here that in some cases it is enough to 
establish domains of existence for homogeneous and inhomogeneous phases. We illustrated 
the method in a model based on a two-dimensional representation of D^h, for which we 
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obtained the following new sufficient conditions: 

For a > the homogeneous state ([|) is stable if and only if \k 2 \ < ^123 + ^/o, < 

«i + ■ 

For a > —j3 and \k 2 \ < K123 + y/@ + cr, < K\+ \ffl + a, irrespective of the sign of a, 
the free energy density is bounded from below. 

Two questions remain open: We don't know the exact domain of stability of the homo- 
geneous phase (|) for a < 0. All we can say is that, as long as eqs. (H) are satisfied, the 
homogeneous phase @ is stable. But phase @ is known to be stable beyond that domain, 
and the transition to the inhomogeneous state is of first order Q. 

The second question concerns the global stability of the free energy beyond the region 
defined by eqs. (fL4l) , which are only sufficient conditions. Consideration of a single plane 
wave 0] imposes some constraints on the coefficients of the free energy, but little is known 
about general configurations. There are, however, two situations where we know that eqs. 
(|I4|) are also necessary conditions for the free energy to be bounded from below: 

For f3 = 0, conditions (0) and (0) are identical, and we can use a scaling argument. 
Beyond the domain defined by eqs. (^4|), the free energy will take a negative value / for a 
modulation of the form ([13]). A similar configuration with V — > nV, A — > nA, and k — > nk 
will satisfy the same boundary conditions for n integer and will have a free energy n 4 f. 
Since n can be arbitrarily large, it is obvious that the free energy is unbounded from below. 

Note that this argument does not apply for f3 7^ 0. But we have studied the system 
numerically and found that, for K23 = ^2 + ^3 = 0, the free energy becomes unbounded from 
below as soon as (|T4] ) is violated. However, our results for k 23 7^ are inconclusive. 
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